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We present the first observation of tiie decay p^K^ , using a data sample of 348 fb~^ 

collected at the T(45') resonance with the BABAR detector. The branching fraction and charge 
asymmetry are measured to be (S.O^j^'g ± 0.6) x 10~* and ( — 12.2 ± 16.6 ± 2.0)%, respectively, where 
the first uncertainty is statistical and the second is systematic. The charge asymmetry is defined by 
^ch = (r^- — ^ B+) I B- + ^B+) with Vg± the decay rate. The significance of the observed 
branching fraction, including systematic uncertainties, is 7.9 standard deviations. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM) of particle physics, the 
weak-current couplings of quarks are described by ele- 
ments of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix [ij. Charmless decays of B mesons provide impor- 
tant information about these couplings. These decays, 
which have branching fractions of the order of 10^^, are 
generally expected to occur via 6 — *■ s or 6 — > d virtual 
loop ("penguin") amplitudes, tree-level b —>■ u decays, or 
a combination of the two. Phenomenological fits to the 
branching fractions and charge asymmetries of charm- 
less B decays can be used to understand the relative im- 
portance of tree and penguin amplitudes and to extract 
measurements of the CKM phase angles. 

We present the first observation of the charmless 6 — > s 
process 5+ p^K^ . Throughout this paper, the charge 
conjugate channel is implied unless otherwise stated. We 
measure the branching fraction and charge asymmetry. 
The latter is defined as Ah = b- b+) I b- +^ b+) 
with r^i the B^ decay rate. A non-zero value of Ah 
implies violation of CP symmetry. Data were collected 
with the BaBAR. detector at the PEP-II asymmetric e+e^ 
collider at the Stanford Linear Accelerator Center. The 
data used in the analysis are based on a sample with 
an integrated luminosity of 348 fb~^, corresponding to 
383±4 million BB pairs recorded at the T{AS) resonance 
[center-of-niass energy (CM) -y/s = 10.58 GeV]. 

The B'^ p^K^ decay is expected to be a pure pen- 
guin decay , making it particularly helpful to separate 
the contributions of tree and penguin amplitudes in other 
channels. Phenomenological studies 0,11, 01 of charm- 
less, strangeness changing {\AS\ — 1) B — > VP decays, 
with V a vector and P a pseudoscalar meson, assume 
that the penguin amplitudes p'y and p'p are related by 
Py = —p'p, where p'y (p'p) is the amplitude for the spec- 
tator quark to appear in the V (P) meson. Measurement 
of the i?+ p'^K^ branching fraction can provide a di- 
rect test of this assumption [5| . Exploiting U -spin sym- 
metry, Soni and Suprun Q recently introduced a tech- 
nique to determine the CKM phase angle 7 with preci- 
sion comparable to the best current measurements, using 
charmless B^ — > M^M" decays, where and are 
charged and neutral mesons. Of the eight M±M° chan- 
nels necessary to apply this technique to B^ —t V^P'^ 
decays, experimental results exist for all but two chan- 
nels: B+ p+K", the topic of this study, and K*+K°. 

Theoretical predictions of the branching fraction for 
_B+ P'^K'^, based on QCD factorization 0,0], heavy 



quark effective theory and flavor SU(3) symme- 
try 0, [l3|, vary from 10~^ to 10~^. The only current 
experimental result is B{B+ p^K") < 4.8 x 10~^ at 
90% confidence level (CL) [ll|. The charge asymmetry 
for this decay is expected to be zero. Any significant de- 
viation from this expectation could provide evidence for 
the creation of non-SM particles produced in the loops. 

The BaBAR detector is described elsewhere In 
brief, charged particle tracks are detected and their mo- 
menta measured by a combination of a five-layer double- 
sided silicon microstrip detector (SVT) and a 40-layer 
drift chamber (DCH), both operating in the 1.5 T mag- 
netic field of a superconducting solenoid. Tracks are iden- 
tified as charged kaons or pions using specific energy loss 
measurements in the SVT and DCH as well as radiation 
angles measured in a ring imaging Cherenkov detector. 
Photons are reconstructed from energy clusters deposited 
in a CsI(Tl) electromagnetic calorimeter. 

Monte Carlo (MC) events are used to determine sig- 
nal and background characteristics, optimize selection 
criteria, and evaluate efficiencies. Samples of e+e^ — > 
r(4S') B°B° and B+B' events, generated by the 
EvtGen [13] event generator, are passed through the 
GEANT-based [IJ] BaBAR detector simulation. The num- 
ber of MC events corresponds to about three times the 
integrated luminosity of the data. We follow a blind pro- 
cedure in which the optimization and systematic study 
of selection criteria, and tests of the fitting procedure, 
described below, are completed before the data are ex- 
amined in the region where the signal is expected. 

A B meson candidate is kinematically characterized by 
the beam-energy substituted mass toes = v's/4 — {p^y 
and the energy difference AE = Ep — •y/s/2, where Ep 
and p*p are the CM energy and 3-momentum of the B 
candidate, respectively. Signal events peak at the nomi- 
nal B mass for toes and at zero for AE. 

We reconstruct i?+ p^K^ candidates through the 
decays K'^ — > tt+tt^ and p~^ 7r"'"7r°, with 

tt" — !• -"y-y. The 7 energy in the laboratory frame is re- 
quired to exceed 30 MeV. The 7r° candidates are required 
to have a mass in the interval [0.115, 0.150] GeV/c^ and 
a laboratory energy larger than 0.2 GeV. The tt" mass 
resolution is about 6 MeV/c^. To improve the resolution 
of Toes and AE, the 7r° candidate's mass is constrained 
to its nominal value. The tt^ candidate is combined with 
an identified charged pion to form a p"^ candidate, which 
is required to have a mass 771.^+^^0 in the interval [0.5, 
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1.0] GeV/c^. The helicity angle Op, defined as the angle 
in the p'^ rest frame between the direction of the boost 
from the _B+ rest frame and the 3-momentum of the 7r+ 
from the p+ decay, is required to satisfy |cos0p| < 0.9, 
since mis-reconstructed tt*^ mesons are concentrated near 
I cos (6'p) I K, 1. Wc form K'^ candidates by combining 
all oppositely charged pairs of tracks, by fitting the two 
tracks to a common vertex, and by requiring the mass to 
lie in the interval [0.490, 0.506] GeV/c^ assuming the two 
tracks to be pions. The mass resolution is about 3 
MeV/ c^. The angle a between the flight direction and 
its momentum vector is required to satisfy cos a > 0.995, 
where the flight direction is the direction between the 
primary and secondary vertices. The X!? candidate is 
combined with the p'^ candidate to form a _B+ candidate 
with a vertex constrained to the beam spot. The 
decay length significance, defined as the ratio of the dis- 
tance between the and 5+ decay vertices and the 
uncertainty on that quantity, is required to be larger 
than 5. The probabilities of the fitted and 
vertices are each required to exceed 0.5%. i?"^ candi- 
dates are required to satisfy 5.25 < tties < 5.29 GeV/c^ 
and |Ai?| < 0.20 GeV. The typical resolution for toes 
(A£;) is approximately 3.0 MeV/c^ (30 MeV). We find 
that 9.8% of the events contain two or more i?+ p^K'^ 
candidates. These are mostly events with more than one 
reconstructed tt". For these events, the candidate with 
the largest B vertex fit probability is retained. 

Backgrounds arise primarily from random combina- 
tions of tracks and clusters in e+e^ ^ QQ [q = u, d, s, c) 
continuum events. To suppress these events, we follow 
a procedure similar to that described in Ref. [13]. We 
use the angle 6t between the thrust axis of the B can- 
didate's decay products and the thrust axis determined 
using the remaining charged tracks and neutral clusters 
in the event, evaluated in the CM frame. The distribu- 
tion of I cos 9t I is nearly uniform for the almost-isotropic 
BB events and sharply peaked near 1 for the jetlike con- 
tinuum events. We require |cos6't| < 0.9. Additional 
use of the event topology is made by employing a Fisher 
discriminant J^, constructed from the angles with respect 
to the beam axis of the B momentum and the B thrust 
axis, and the energy fiow around the B thrust axis 15 1. 

Potential backgrounds from B^B~ and B'^B'^ events 
arise from B Dtt, K*{892)tt and K^{U30)tt decays 
that have the same 7r+7r°ifg final state and similar peak- 
ing structure in toes and AE as the signal events. The 
selection requirement applied to m^+^o, given above, re- 
jects most of these backgrounds. To further reduce the 
BB background, we apply a veto (1.78 < to^o^o < 
1.94) GeV/c^, a D+ veto (1.83 < to^o^+ < 1.91) 
GeV/c2, a K*° veto (0.8 < m^o^o < 1.0) GeV/c^, a 
K*+ veto (0.8 < TO^o^+ < 1.0)''GeV/c2, a /s:*"(1430) 
veto (1.3 < TO^o^o < 1.6) GeV/c^, and a i4:*+(1430) 
veto (1.3 < mj^o^+ < 1.6) GeV/c^, where a veto indi- 



cates that an event is rejected if the two-particle invari- 
ant mass lies in the specified mass window. These veto 
criteria are determined as follows: 4 standard deviations 
of the experimental resolution around the mass peaks for 
D'^ and D", and 2 (1) resonance widths PJI] around the 
mass peak for K* [if* (1430)]. The decay B+ K°Tr+ 
{B^ — * -^s""") can contribute to background when the 
decay products are combined with a low momentum tt'^ 
(7r+). We therefore require the K^-k^ and K^-k^ invariant 
masses to be less than 5.2 GeV/c^. 

These criteria reject more than 99% of the BB back- 
ground channels discussed above and about 30% of the 
signal. The remaining BB background is combinatoric 
and does not peak in AiJ and toes- 

From MC simulation, the signal efficiency is deter- 
mined to be (14.78 ± 0.10)% where the uncertainty is 
statistical. This efficiency has been corrected to account 
for small differences in neutral particle reconstruction ef- 
ficiencies between the data and MC, and for differences 
in the identification efficiency of the 7r+ used to recon- 
struct the p^ . As an example, the latter correction is 
determined using a D*'^ Z?°7r+ data control sample 
with Z?" K~TT~^. The efficiency corrections are 97% 
for the 7r° reconstruction and greater than 99% for the 
K'^, reconstruction and 7r+ identification. 

The number of signal events (the signal yield) and 
charge asymmetry are determined from an extended un- 
binned maximum likelihood (ML) fit with the following 
variables: toesj Ai?, JF, to^+^o, cosOp, and the B flight 
time significance, with the latter variable defined as the 
proper time difference At between the produced B and 
B candidates divided by its uncertainty (ta* [17|. The 
B vertex is determined by fitting all tracks except the 
daughters of the signal B candidate to a common vertex, 
employing constraints from the beam spot. The likeli- 
hood function has the form 



^-^exp(- 



11% 



N 

n 

i=l 



^71j-Pj(x,) 

J = i- 



(1) 



where N is the total number of input events, Uj is the 
fitted yield of component j (signal, continuum, and BB 
background), and Vj{'x.i) is the corresponding overall 
probability density function (PDF), given by 



xVj{cosep)rj{At/aAt)- 



(2) 



The signal and BB background PDFs are determined 
from MC simulation. The continuum background PDF 
is obtained from sideband data (0.1 < \AE\ < 0.2 GcV 
for Toes and 5.25 < toes < 5.27 GeV/c^ for other vari- 
ables). For toeS) the PDFs of the signal and contin- 
uum are parameterized by a Crystal Ball [l8| and an 
ARGUS function [l^, respectively. A relativistic Brcit- 
Wigner function with a p-wave Blatt-Weisskopf form fac- 
tor [20| is used to model the signal TOtt+ttO distribution. 
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For the background components, the 7r+ tt*^ mass is mod- 
eled by a combination of a polynomial and the signal 
function. Slowly varying distributions {AE for the con- 
tinuum background, and cos^p) are modeled by polyno- 
mials. The remaining variables are parameterized with 
either a Gaussian, the sum of two or three Gaussians, or 
an asymmetric Gaussian. Dips occur near | cos9p\ = 0.81 
because of the resonance vetoes. We describe these dips 
by two Gaussian shapes. We use a large data control 
sample of B+ D°Tr+ {D" K^tt°) events to verify 
the simulated resolutions and peak positions of the rriEs, 
AE, T, and At /a At signal PDFs. 

Eq. ^ is based on the assumption that the variables 
in the PDFs are uncorrelated. To evaluate possible bias 
in the signal yield that might arise from residual correla- 
tions, we construct an ensemble of 600 simulated exper- 
iments. Each experiment contains the expected number 
of signal, continuum, and BB background events. The 
continuum events are randomly drawn from the PDFs 
while the signal and BB background events are randomly 
drawn from the MC samples. The bias is defined as 
the difference between the mean signal yield, determined 
from fits to the simulated experiments, and the number 
of signal events included in the samples. The bias in the 
signal yield is determined to be 4.8 ± 1.2 events, where 
the uncertainty is statistical. 



TABLE I: Summary of results. The uncertainties on the event 
yields, fit bias, and efficiencies are statistical only. 



Parameter 


Value 


Events in fit 


41150 


Signal yield (events) 


1581^^ 


Continuum background yield (events) 


40,321^^^° 


BB background yield (events) 




Fit bias (events) 


4.8±1.2 


Detection efficiency (%) 


14.78±0.10 


Daughter branching fractions Y[ Si (%) 


34.18±0.03 


Statistical significance (a) 


8.2 


Significance with systematics (a) 


7.9 


Branching fraction B (xlO"'') 


8.01^ * ±0.6 


Charge asymmetry (%) 


-12.2 ±16.6 ±2.0 



Table U lists the results of the fit to the data. The fit 
yields a simultaneous determination of the fraction 
(/_) of 5+ (-8^) events relative to the total number of 
signal events, with the constraint /+ ± /- = 1. The 
charge asymmetry is determined from = 2/_ — 1. 
The statistical uncertainty of the signal yield is given 
by the change in the central value when the quantity 
— 21n£ increases by one unit from its minimum value. 
The statistical significance is given by the square root of 
the difference between the value of —2 ln£ for zero signal 
events and the value at its minimum. The correspond- 
ing significance including systematic uncertainties (dis- 
cussed below) is determined by convolution of the likeli- 
hood function with a Gaussian distribution whose stan- 



dard deviation equals the total systematic uncertainty. 
Figure [T] shows projections of the fitted variables. To en- 
hance the visibility of the signal, events are required to 
satisfy C^{S)/[C^{S) + Ci{B)] > 0.9 (this retains J70.0%, 
1.4%, and 14.5% of the signal, continuum, and BB back- 
ground events, respectively), where £i{S) is the likeli- 
hood function for signal events excluding the PDF of the 
plotted variable i and Ci{B) is the corresponding sum for 
all background components. 




0.6 0.8 1 ^05 ' " 0^5 

m,^,,. (GeV/c^) cosBp 




FIG. 1: Distributions of (a) AE, (b) tties, (c) m^+^o, 
(d) cosOp, (e) T , and (f) At/aAt- To improve visibility, a se- 
lection requirement on the likelihood ratio that retains 70.0% 
of the signal events has been applied. The points with uncer- 
tainties are the data. The curves are projections of the ML 
fit. The dashed curves show the sum of the continuum and 
BB background components. The dot-dashed curves show 
the signal component. The solid curves show the sum of the 
signal and background components. 

We calculate the branching fraction by subtracting the 
fit bias from the measured signal yield and dividing the 
result by the overall efficiency and the number of pro- 
duced BB pairs N^^. The overall efficiency is the prod- 
uct of the detection efficiency and the daughter branch- 
ing fractions [lB| (see Table |T|. We assume equal decay 
rates of the r(45) to B+B' and B°B°. The branch- 
ing fraction and charge asymmetry are determined to be 
iS.Otli ± 0.6) X 10-6 and (-12.2 ± 16.6 ± 2.0)%, respec- 
tively, where the first uncertainty is statistical and the 
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second is systematic. We determine —0.40 < Ach < 0.15 
at 90% CL, including systematic uncertainties. 

The principal sources of systematic uncertainty are as 
follows. The uncertainty in the tt^ reconstruction effi- 
ciency is 3%. The uncertainty related to the signal PDFs, 
assessed by varying the fitted PDF parameters within 
their uncertainties as determined from the 5+ Z?°7r+ 
data control sample, is 3.2 events. An uncertainty in the 
fit bias (2.7 events) is defined by the quadratic sum of half 
the bias itself and the statistical uncertainty of the bias. 
To evaluate the effect of a possible non-resonant compo- 
nent, we generate a Monte Carlo sample using a 3-body 
Dalitz amplitude event generator for — > 7r+7r°iCg, in- 
cluding the if * (892), if 5 (1430), p resonances and a non- 
resonant amplitude. Re-performing the ML fit with sig- 
nal PDFs determined from this sample results in a 3.5% 
increase in the signal yield, which we take to be the sys- 
tematic uncertainty. Variations of all resonance vetoes 
yield an uncertainty of 3.1%. When the requirement on 
cos^p is varied, the results change by 2.0%. Other prin- 
cipal sources of uncertainty are those from the track re- 
construction efficiency (1.6%), the BB background PDFs 
(2.0 events), N^g (1.1%), and variation of the selection 
criteria on | cos^tI (1.0%). We add all terms in quadra- 
ture to obtain the total systematic uncertainty. 

In summary, we present the first observation of the 
pure penguin b ^ s decay process i?+ p'^K^. The sig- 
nificance of the measured branching fraction is 7.9 stan- 
dard deviations. Using the assumption p'y — —p'p @, 
the i?+ p^K^ branching fraction is predicted to lie 
between about 9 and 13xl0~^, consistent with our mea- 
surement within the uncertainties. The measured charge 
asymmetry is consistent with the SM expectation of zero. 
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